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Abstract 



The integral representations for the eigenfunctions of N particle quantum open and 
periodic Toda chains are constructed in the framework of Quantum Inverse Scattering 
Q-i! Method (QISM). Both periodic and open iV -particle solutions have essentially the 

same structure being written as a generalized Fourier transform over the eigenfunctions 
of the N — 1 particle open Toda chain with the kernels satisfying to the Baxter 
equations of the second and first order respectively. In the latter case this leads to 
recurrent relations which result to representation of the Mellin-Barnes type for solutions 
of an open chain. As byproduct, we obtain the Gindikin-Karpelevich formula for the 
Harish-Chandra function in the case of GL(N, K) group. 
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1 Introduction 

This report is devoted to well-known quantum mechanical problem to find the simultaneous 
eigenfunction of commuting set of Hamiltonians for the periodic Toda chain. The first 
important step in this direction has been done by Gutzwiller |] who solved the problem 
for the particular cases N = 2,3 and 4 particles and found such important phenomena 
as quantization of spectrum and separation of multidimensional Baxter equation into the 
product of one dimension ones. In fact, he performed the quantization of periodic Toda chain 
in terms of separated variables introduced by Flaschka and McLaughlin ||. Next important 
step was taken by Sklyanin || who constructed R -matrix formalism for both classical and 
quantum cases Toda chains and introduced the algebraic method of separation variables for 
an arbitrary number of particles. His approach drastically simplifies the derivation of the 
Baxter equation and works for wide spectrum of integrable models |4| . 

Our method to solve the spectral problem consists of analytical re-interpretation of 
Sklyanin's algebraic ideas which allows to find the integral representation for the eigen- 
functions of the periodic Toda chain as a kind of generalized Fourier transform with the 
eigenfunctions for the open Toda chain ||. In turn, this method can be treated as a natural 
generalization of original Gutzwiller's approach. The explicit solution for the eigenfunctions 
of the open Toda chain plays a key role in this construction. 

It has been discovered by Kostant |J that the commuting set of Hamiltonians of an open 
Toda chain coincides with the Whittaker model of the center of universal enveloping algebra. 
Hence, the Whittaker functions are in fact the eigenfunctions for the open Toda chain. In the 
usual group-theoretical way the Whittaker function is defined as a matrix element between 
compact and Whittaker vectors 0-@ in the principal series representation. There are many 
obstacles to generalize this approach to other quantum models or to loop groups. 

The present approach to construct the eigenfunctions for both periodic and open chains 
is rather different || [l(| : it bases on the Quantum Inverse Scattering Method for the peri- 
odic Toda chain ||. One of the interesting results of analytical calculations in the R -matrix 
framework is the revealing of recurrent relation between N and N — 1 particles eigenfunc- 
tions for the open Toda chain (in fact the idea to use a recurrent relation was pointed out 



by Sklyanin in 11]; our recurrent relations is an explicit realization of such an idea). This 
naturally leads to new integral representation for the Weyl invariant Whittaker functions to 
compare with classical results 0-[[J. This representation is quite explicit and very useful to 



investigate the different asymptotics. In particular, the Gindikin-Karpelevich formula p2 



for the Harish-Chandra function |T3| can be obtained in a very simple way for particular 
case of GL(N, M) group. The eigenfunction for the periodic Toda chain are constructed in a 
rather explicit form and have essentially the same form as the recurrent relation mentioned 
above. The integral formula for eigenfunctions can be considered as a representation of the 
Whittaker functions for GL(N) group at the critical level. 

The present approach can be generalized to other quantum integrable models. For exam- 



ple, the eigenfunctions for the relativistic Toda chain are calculated in |l y using the same 
QISM ideology. 



2 Quantum Toda chain: description of the model 

2.1 Periodic spectral problem 

The quantum iV -periodic Toda chain is a mult i- dimensional eigenvalue problem with N mu- 
tually commuting Hamiltonians Hk(xi,pi; . . . ; x N ,p N ) , (k = 1, . . . , N), where the simplest 
Hamiltonians have the form 

N 

Hi = "^Pk 
fc=i 

N 

H 2 = j2p*p™-J2 eXk ~ Xk+1 (2,1) 

k<m fc=l 

H 3 = ^ PkPmPn + ■■; 
k<m<n 

(x N+ i = xi) etc. and the phase variables Xk,Pk satisfy the standard commutation relations 
[xfe,p m ] = ih5k m ■ The main goal is to find the solution to the eigenvalue problem 

H k % = E k \ k = l,...,N (2.2) 

with fast decreasing wave function \l/ . To be more precise, let us note that, due to trans- 
lation invariance, the solution to ( |2.2j ) has the following structure: 

(i N } 
^ E (xt,...,x N ) = $ E {x 1 -x 2 , ...,x N -i-x N ) exp< -Ei^x k > (2.3) 



fc=i 



One needs to find the solution to (|4.11|) such that *$? E G L 2 (W N l ) . In equivalent terms, we 
impose the requirement 

J f{E x )^ E {x x , . . . , x N )dE 1 G L 2 (R N ) (2.4) 

for any smooth function f(y) , (t/6R) with finite support. 

2.2 GL(N - 1,R) spectral problem 



It turns out that solution to (|2.2j ), (|2.4j ) can be effectively written in terms of the wave 
functions corresponding to open N — 1 -particle Toda chain (quantum GL(N— 1,R) chain). 
The Hamiltonians of the latter system can be (formally) derived from (|2.1| ) by cancelling 
out all the operators containing p N and x N thus obtaining exactly N — 1 commuting 
Hamiltonians h k (xi,pi] . . . ;x N -i,p N ^i) (k = 1,. . . , N—l). Let 7 = (71, . . . ,7^-1) € R 7 ^ 1 , 
x — (xi, . . . , x N -\) G R^ -1 . We consider GL(N — 1, R) spectral problem 

M<r(*) = ff*(7)^(*) fc = l,...,JV-l (2.5) 



where &k{l) are elementary symmetric functions. 

Obviously, in the asymptotic region Xk+i 3> Xk, (k = 1, . . . , N—2) all potentials vanish and 



the solution to ( |2.5|) is a superposition of plane waves. The problem is to find a solution to 



75]) satisfying the following properties: 
(i) The solution vanishes very rapidly 

^{x) ~ exp { - J* e (* h -**+i)/2j Xk - Xk+1 _> oo (2.6) 

(ii) The function -0 7 is Weyl-invariant, i.e. it is symmetric under any permutation 

i>... 7j ...y k ... = ^... 7fc ...7,... ( 2 - 7 ) 

(iii) ipy can be analytically continued to an entire function of 7 G C^ -1 and the following 
asymptotics hold: 

V7 ~ |7,l^exp{-^(iV-2)| 7j |} (2.8) 

as |Re7j| —>■ 00 in the finite strip of complex plane. 
The properties (i)-(iii) define a unique solution to the spectral problem (|2.2|). 



3 Main results 

Theorem 3.1 The following statements hold /|3|, [IT]/.- 

(i) Let a set \\~Yjk\\ be the lower triangular (JV— 1) x (JV— 1) matrix. The solution to 
the spectral problem flg.^ - pTffi ) can be written in the form of multiple Mellin-Barnes 



integrals Q: 



ih 



Uk=i 2 k\ J t\ n r( 7 ^~ 7 " fc )r( 7 " fc ~ 7 ^) 



i,fe=i 



X 

(3.1) 



JV-l "1 N-2 



I . iV— 1 I JV— z 

eXP 1 fc 5Z Xn ( 7nfc ~ ln-l,k) } IJ ^7j 



n,k=l ) j,k=l 

where the integral should be understand as follows: first we integrate on 711 over the 
line Ini7n > max{Im7 2 i, 1UL722} ; then we integrate on the set (721, 722) over the lines 
1H172J > max m {Im7 3m } and so on. The last integrations should be performed on the 
set of variables (7^-2,1 • • • , 7^-2,^-2) over the lines lmj N _ 2t k > max m {Im7 JV _ 1]m } . 



1 We identify the set 7 with the last row (7^—1,1, ■ • ■ ,7^-1,^-1) 



(ii) In the region Xk <C Xk+i{k = 1, . . . , N — 1) the solution has the following asymptotics: 
^{x) = J2 <P{si)e^ s " f ' x) + 0(max|e Xfc - a;fe + 1 } iV_1 ) ( 32 ) 



sew 



where (.,.) is a scalar product in M. N x and the summation is performed over the 
permutation group; 0(7) is (renormalized) Harish- Chandra function 

m = n - *,,>/» nr(^) (3. 3) 



j<k 

JV-1 

where (7, p) = | J] (^ _ 2771)7*. ■ 

m=l 

(mj T/ie functions t \3. 1\ ) have the scalar product 

s&W 

and obey the completeness condition 

p(j)^(x)^(y)dr( = 8(x - y) ,^ 

where 



J ^,(x)^(x)dx = i£M J2 5 ('^ - V) (7,Y 6 k^ 1 ) 



(3.4) 



is the Sklyanin measure fi^J. 



l-N .- - . -2 



7j~7fc 
ih 



(3.6) 



The eigenf unctions for the periodic chain are constructed as a kind of Fourier transform with 
the function ( |3.1|) . Let 



TV 



t N (\;E) = J2(-V k * N ~ k E k (3.7) 

and Bj denotes j -th basis vector in 



fc=0 
)JV-1 



Theorem 3.2 T7ie solution to the spectral problem ( \2.2[ ), ( \2.4\ ) can be represented as the 
integral over real variables 7 = (71, . . . , 7jv_i) in the following form: 

® E (x,x N ) = — / i2(-r)C(T,E)^^^ El (x,x N )d-f ^^ 

where 



(i) The function ^~f,Ei (x, x N ) is defined in terms of solution ( W. 1\ ) to the GL(N — 1,M.) 
spectral problem: 

JV-l 

Vy tEl {x,Xir) = if> y (x) exp|-^i-^7 m ja; J v| (3.9) 

m=l 

(ii) The function C(*y; E) is the solution of multi-dimensional Baxter equations 

t N ( lr , E)C(r, E) = i N C(>y + ihef, E) + i~ N C(-y - ihef, E) (3.10) 

which is symmetric entire function in 7 -variables with the asymptotics 

C(7;S)~| 7fe |^/ 2 exp{-^M} (3.11) 

as Re7fc — > ±00 in the strip |Im7fc| < H 



The above restrictions imposed on solution to ( |3.10|) are reformulation of the quantization 



condition ( |2.4p on the level of 7 -representation. To obtain the explicit integral form for the 
eigenfunctions, we use the solution to fl3.10|) , ( [3.1 H ) i n the Pasquier-Gaudin form |L5] (see 
sect. 7 below) 

C ^> E > = 11 — ; r- (3.12) 

6 k (EV 



J. iV f 

=1 nsmhf ( 7i 

fe=i v 



where the entire functions c±('y) are two Gutzwiller's solutions [|TJ of the one- dimensional 
Baxter equation 

£(7; E)c(r, E) = C N c(-y + ih; E) + 1*0(7 - ih; E) (3.13) 

and the parameters £(E), 5 = (61(E), . . . , S N (E) ) satisfy the Gutzwiller conditions (the 
energy quantization) Jl|, jl5f(see sect. 7 below). Then the multiple integral (|3.8| ) can be 
explicitly evaluated. Let y = (yi, . . . ,y N ) G M. N be an arbitrary vector. We denote 
y^ = (yi, . . . , y s -i, y s +i, • • • , Vn) the corresponding vector in M^ -1 . 



Theorem 3.3 Assuming that 6j(E) 7^ 6k(E) , the solution (3.8) can be written (up to an 
inessential numerical factor) in the equivalent form 

N 

%(x,x N ) =£(-1)""' YI ^(S is) +ihn^)C + (S^+ihn^)^ si s )+ihn(s)>Ei (x,x N ) 

8=1 n^ez 1 ^- 1 

(3.14) 

where 

JV-1 

C+(<y) = J] c+faE) (3.15) 

i=i 
and A(7) = YI (lj ~ Ik) is the Vandermonde determinant. 

3>k 



Remark 3.1 For N = 2, 3 and N = 4 formula ( 3.14 ) reproduces the results obtained by 
Gutzwiller Fi|/. 



4 R -matrix approach 

The Toda chain can be nicely described using the R -matrix approach |J. It is well known 

that the Lax operator 

f X-p n e~ x " \ 

Ln(X) = . (4-1) 

V -e Xn / 

satisfies the commutation relations 

R(X - fi)(L n (X)) ® /)(/ ® L n (^)) = (/ ® L n (/i))(L n (A) ® I)R(X - pi) (4.2) 

where 

J R(A) = /®/ + ^P (4.3) 

A 

is a rational R -matrix. The monodromy matrix 

r " (A)£MA) - Ll(A)S ( C „(A)^(A)J (44> 

satisfies the analogous equation 

R(X - fi)(T(\) ® /)(/ <g> TO*)) = (/ ® T(n))(T(X) ® I)R(X - /i) (4.5) 

In particular, the following commutation relations hold: 

[A W (A), ^(/i)] = [C N (X), CM] = (4.6) 

(X- fx + ih)A N (fx)C N (X) = (X- fi)C N {X)A N (fx)+ihA N (X)C N (fi) (4.7) 

(A - fi + ih)D N {X)CM = (A - h)C n (h) d n(X) + iflD N {pi)C N {X) (4.8) 

From ( |4.5|) it can be easily shown that the trace of the monodromy matrix 

t N (X) = A N (X) + D N (X) (4.9) 

satisfies the commutation relations [t(A),t(/x)] = and is a generating function for the 
Hamiltonians of the periodic Toda chain: 

N 

?*(A) = £(-l)*A*-*ff fc (4.10) 

fc=0 

We reformulate the spectral equations ( |2~2| ) as follows: 
where 

JV 

t N (X;E) = J2(-±) k * N - k E k (4.12) 



fc=0 



6 



On the other hand, it can be easily shown that the operator 



7V-1 



A W _ X (A) = ^(-l)^-*- 1 /^!,^;...;^,^) (4.13) 



fc=0 



is nothing but the generating function for the Hamiltonians hk of GL(N-l) Toda chain. 
Therefore, the GL(N — 1,IR) spectral equations can be written in the form 



JV-l 



A N „ x (\)i; y (x) = \\{\- lm )^{x) (4.14) 



m=l 



Using the obvious relation 

C N {\) = -e XN A N ^(X) (4.15) 

one obtains, as a trivial corollary of (|4.14j ), 

C N (ij)^(x) = V 7j e 7 (4.16) 



Remark 4.1 Equations (\4-H\ ) are an analytical analog of the notion of "operator zeros' 
introduced by Sklyanin j^j. 



5 Eigenfunctions for the open Toda chain 



Suppose that solution to (|4.14j ) satisfying ( p.6|) -( p78|) is given. Using the commutation rela- 



tions (|4.7j), ( |4.8D together with ( |4.1t>| ), it is easy to show that the following relations hold 

^iv(7i)^ = i~ N e~ XN ip~ f - i M J (5.1a) 

D N {lM~i = i N e XN ip-y +i ^ (5.1b) 

(j = 1, . . . ,N — 1) where e.j is j-th basis vector in R^ -1 . Note that ( |5.1b| ) is a corollary 
of ( |5.1a| ) since the quantum determinant of the monodromy matrix ( |4.4|) is unity. 



Let us introduce the key object - the auxiliary function 



2V-1 



$? y>e (xi,...,x N ) = ^ 7 (a;)exp|-^e -^7 m jx JV | (5.2) 



m=l 



where e is an arbitrary parameter. From ( [4.14| ), ([4.15 ) and (|5.1| ) it is readily seen that this 
function satisfies to equations 

AT-l 

C N (A)^, 6 = - e XN H (A - Ti ) * T , e (5.3a) 



N-l N-l N-l 



A„(\)%, e = (A-e + ^ 7m ) H(\-jj) ^e + i^Y,^-^ \\ —^ (5.3b) 

A-7n 
j=l mjtj 

In particular, 



m=l J=l J = l m 7v 

iV-1 

^(A)* xe = i" £ ^+iBe i)e J] ^3^ ( 5 - 3c ) 

i=l m# 7 -* 7 " 



^(7,)*^ = ^^,e + r"* v (5.4) 

The problem is to find the corresponding solution for GL(N, R) Toda chain using the above 
information, i.e. in terms of the function ty^ tE (x) construct the Weyl invariant function 
V'Ai,...,Ajv( x i> ■ ■ ■ > x iv) satisfying to equations 

AT 

A N (X)ipxi,...,x N = \\(X - A fe ) ^Ai,...,Aiv ( 5 - 5a ) 

Jfe=l 

Av +1 (A n )^ Al ,..., A]V = r^" 1 e XN ^x u „.,x n -in,...,x N (n = 1, . . . , N) (5.5b) 



and obeying the similar to (|2.6|)-(|2lf) conditions 



Lemma 5.1 /|7^/ Let tyy je (x,x N ) be the auxiliary function (\5.3Q . Let A = (Ai, . . . , X N ) € C N 
be the set of indeterminates. Let 

(27rh) 1 - N yj j /7j-7fc\ r flk-lj 



*)^Srn^)rm (5. 6 ) 



j<k 

N-l N 



g( 71 ,...,7.-iiAi,...,A,)=nn^'^ r ( 2 ^) ( 5 - 7 ) 

j=l k=l 

Then the Weyl invariant solution to the spectral problem ( 5.5dj) -( '5.5l ) with the properties 
discussed above is given by recurrent formula 

ipx 1 ,...,x N {xi,...,x N ) = / fjt(i)Q('y;X)^ ; x 1 +...+x tr (xi 1 ... 1 x If )d'y ( 58 ) 

c 

where the integration is performed along the horizontal lines with Im^j > max/,, {Im A&} . 

Proof. One needs to calculate the action of the operators A N (X) and 

A N+1 (X) = (A - p N+1 )A N (X) + e~ x ^C N (X) (5.9) 



on the function (|5.8|) using the formulae fl5.3b|) and ( 5.5a|) , ( |5.3a| ). The shifted contours can 



be deformed to original ones using the facts that integrand in ( |5.8| ) is an entire function fast 
decreasing in any finite horizontal strip of complex plane as | Re 7^ | — >• 00 . The last step is 
to use the difference equations for the parts of the integrand with respect to the shifts ±ih 

8 



of parameters 7™ and \k ■ 



Proof of Theorem |3.1| . The proof of (|3.1| ) is straightforward resolution of the recurrent 



relations (|5.8|) starting with trivial eigenfunction ip 1 (xi) = exp{|7nXi}. Obviously, the 
function ( [3.1|) is symmetric under the permutation of parameters 7 . The asymptotics fl2.6|) 
can be proved using the steepest descent method. Using the Stirling formula for the T- 
functions as 7^-1, k = Ik —> ±00 , it is easy to see that the asymptotics (|2.8| ) hold. Hence, 
( |3.1| ) is an appropriate solution to the spectral problem. 

Further, the formula ( |3.2|) can be proved as follows. The integrand in ( |5.8|) decreases expo- 
nentially as |tj| — > 00 , (j — 1, . . . , N— 1) in the lower half-plane and, as consequence, the 
integrals over large semi-circles in the lower half-plane vanish. Using the Cauchy formula 
to calculate the integral (|5.S| ) in the asymptotic region x fc+1 3> x^, {k = 1, . . .,N — 1) , it 
is easy to see that the asymptotics of the function ip^ are determined precisely in terms of 
the corresponding Harish-Chandra function (p.3|). 



The scalar product ( |3.4| ) is the consequence of the Plancherel formula proved in |16| for the 
SL(N,M.) case. The formula ( p75|) can be proved by induction. ■ 



Remark 5.1 In ||/ (eqs. (4-7), (4-18)) the eigenf unctions for an open Toda chain have been 
constructed in terms of Whittaker function which has a standard integral representation 
corresponding to the Iwasawa decomposition of semi-simple group (see for example j^). Our 
expression t \3. 1\ ), being obtained in the framework of Quantum Inverse Scattering Method, 
seems quite different. Nevertheless, both representations do define the same function (this 
can be shown by comparing the corresponding asymptotics and analytical properties). Hence, 
one can consider the representation ( $.1\ ) as a new one for Whittaker function /Q/-/0/- 



6 Periodic chain: 7 -representation, eigenfunctions, 
and Plancherel formula 



Let \Pg, (x, x N ) be the fast decreasing solution of the problem (|4.11|) . We define the function 
C(7; E) by the generalized Fourier transform: 



5(E 1 — e) (7(7; E) = / & E {x, x N ) ^-y je (x , x)dxdx t 



(6.1) 



BiV-l 



Lemma 6.1 The function C(~f) possesses the following properties: 

(i) It is a symmetric function with respect to 7 -variables, 
(ii) It is an entire function of 7 G C w_1 . 
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(iii) The function C('y) obeys the asymptotics 

C(r,E)~\ lk \- N />e W {-^M} (6.2) 

as Re7fc — > ±00 m the strip | Im ^ I < ^■ 
(raj 27ie function C("y) satisfies the multi-dimensional Baxter equation 

£( 7j ; E)C(r, E) = t N C(-y + ihef, E) + i~ N C(<y - ihef, E) (6.3) 

where £(7; E) is defined by (j^TT^). 



Proof. The symmetry of the function C(-y) is obvious. 
We present here only a sketch of the proof of the statements (ii) and (iii). 
The statement (ii) follows from the assertion that the auxiliary function \I/ 7j e is an entire 
one while the solution of the periodic chain vanishes very rapidly as \x k — Xk+i\ ~~ *■ °° • 0- 
(iii) The asymptotics ( |6.2|) is a combination of two factors. The first one comes from the 
asymptotics (|2.8|) while the additional factor ~ |7fc| _1 exp{— n\^ k \/K\ results from the sta- 
tionary phase method while calculating the multiple integral including the function ( |3.1| ). 
The calculation is based heavily upon the exact asymptotics of the function ^/ E (x,x N ) as 
\x k - x k+ i\ —* 00. 

The proof of (iv) is simple. Using the definition ( [4.11|) and integrating by parts (evidently, 
boundary terms vanish), one obtains 

5(E 1 -e)t(-/ j ;E)C(j)= / |f(7 i )^(a;,x JV )| ^ >e (x,x N )dxdx N = 

^ (6-4) 

^ E (x, x N ) t(~/j)i£?y te (x, x N )dxdx N 

R N-1 

Taking into account the relation fl5.4| ), the Baxter equation ( |6.3| ) follows from definition 

TTD. ■ 



Now we prove Theorem FO. Using the completeness condition 



/ /j,('y)y~ fte (x,x N )y~ f , e (y,y N )d'yde = 2irH8(x - y)S(x N - y> 



(6.5) 



which is a corollary of (|3.5|) , the inversion of the formula ( |6.1|) results to expression 
* E (cc,x JV ) = — / n(j)C(~f,E)^^ jEl (x,x N )d>y 



(6.6) 



2 Actually, the boundary conditions have the same importance here as the requirement of compact support 
in the theory of the analytic continuation for the usual Fourier transform. 
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The integral ( |6.6| ) is correctly defined. Indeed, the measure ( |3.6| ) is an entire function. 



Therefore, there are no poles in the integrand. Moreover, 

M7)~l7 fc r- 2 exp{^(iV-2)|7 fc |} (6.7) 



as |7fc| — > oo. Taking into account the asymptotics ( |2.8|) and (|6.2| ) one concludes that the 
integrand has the behavior ~ |7fc| _1 exp{ — ir\^k\/K\ as \^k\ — > oo . Therefore, the integral 
( |6.5| ) is convergent. 
One can directly prove the spectral problem ( f4.11| ) calculating the action of the operator 



tjv(A) = A N (\) + -Djv(A) on the right hand side of ( |6.6| ) with the help of the formulae ( |5.3b| ) 
( |5.3c|) . The calculation is performed similarly to those of Lemma |5.1|, using the analytical 



properties of the integrand and the Baxter equation ( |B~3|) (see || for details). 



The last step is to prove that the function (|6.6j ) satisfies to integrability requirement ( |2.4| ). 
Using the scalar product 

^^ je/ (x,x N )^^ e {x,x N )dxdx N = (2nh)— — — 5(e-e) J^ 5(«7 ~ V) (6 8) 

one can write the Plancherel formula 

2irh / ^ E ,(x, x N )^ E (x, x N )dxdx N = 8(E 1 -E[) / ^(7) C(t, E')C(t, E)d-y 



(6.9) 



pJV 



The integral in the r.h.s. of (|6.9|) is absolutely convergent due to asymptotics ( |6.2j ) and 
( |6.7|) . Hence, the norm H^H is finite modulo GL{1) 5 -function S(E 1 — E[) (see the 
corresponding factor in ( |2.3| ) which leads to this function) and requirement ( |2.4| ) is fulfilled. 
Hence, Theorem p.2| is proved. ■ 



7 Solution of Baxter equation 

It is well known pi |T3[ (see also M for details) that the solution to the Baxter equation 
( |3.10| ) with the asymptotics ( |3 . 1 1| ) can be written in the following separated form: 



C( T ,E) = li C ^ E) -^-^ E 



N TT (7- 1 ) 

1 nsinh|( Ti -4) 



fc=i 
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where £ and 5k are an arbitrary constants and the entire functions c±( r y) are defined it 
terms of N x N determinants: 



c+(7) 



1 



N _£7 

Y[h T(l-£( 7 -A fc )) 

fc=i 



1 


t(-f + ih) 







1 


1 


1 


n 


t(j+2ih) 


i(7+2ift) 







1 


1 


1 



i(7+3ift) 



£(7+3ift) 



(7.2a) 



c_(7) 



1 



JV £7 

nnsr(l + i( 7 -A & )) 
fe=i 



t(7-3ifi) 


1 


t{~j-Zih) 





n 


1 


1 


1 




t(j-2ih) 


t('y-2ih) 







1 


1 



(7.2b) 



t(-y-ih) 

and Afc = \k(E) are the roots of the polynomial £(7) = ^(7; -E) . 

On the other hand, the solution ( [7. 1|) is not an entire function in general since the denomi- 
nator in ( [r.l|) has an infinite number of poles at 7 = 5k + z/in& , n^ £ Z, A; = 1, . . . , N . The 
poles are cancelled only if the following conditions hold: 



c+(S k + ihn k ) = £c_(4 + ihn k ) 
In turn, this means that the Wronskian 



(7.3) 



(7.4) 



vanishes at 7 = 5k + ihn k . The Wronskian is ^-periodic function and possesses exactly 
N real roots 5 k (E) |J. Therefore, the solution ( |7. 1| ) has no poles if one takes 5k = 5 k (E) 
provided that the constant £ is chosen in such a way that 



£ 



c+(7) 



C-(t) 



jfc 



iV 



7=4 (£) 



(7.5) 



Hence, one arrives at the following 
Lemma 7.1 / fJ^/ 27ie function 



N-l 



c(r, E) = J] 



1 nsmh|( 7j -^(E) 
fc=i ^ 



(7.6) 



where 5 k (E) are real zeros of the Wronskian (\7.4 ) and the constant £ is chosen according 



to ([7.3j) , satisfies to conditions of Lemma \6.1. 
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The quantization conditions 

c +0i) c+{6 N ) 

c_(*!) •" c46 N ) { '- l) 

determine the energy spectrum of the problem. They have been obtained for the first time 
by Gutzwiller [Q using quite different method. 



To prove Theorem |3l|, one should substitute the solution (|7.6| ) into the integral formula 



1|) and calculate the residues coming from individual terms 

c±{lj\E) 



N , , 

nsinh|( 7j -4(^)) 



(7- 



The result is exactly the sum over all possible poles of expressions ( [7.81) and essentially 
coincides with ( |3.14p (see careful analysis in [|5J). ■ 
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